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ABSTRACT 


A  numerical-graphical  method  is  described  to  extract  weak  mesoscale 
signatures  from  Seasat  altimeter  data.  Results  are  presented  in  a  form 
suitable  for  overlay  on  a  Mercator  grid.  The  signatures  presented  are  thought 
to  be  due  to  tectonic  anomalies. 
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INTRODUCTION 

Low  pass  filtering  of  satellite  altimeter  data  to  estimate  sea  surface  topography  and  slope  has  been 
discussed  by  various  authors  [  1 ,2,31 .  These  methods  tend  to  smooth  out  some  physically  meaningful 
signatures  embedded  in  the  data,  especially  those  of  small  scale.  In  order  to  bring  out  these  signatures 
that  are  not  far  above  the  noise  level,  it  is  useful  to  pass  the  data  through  a  high  pass  filter.  The  digital 
filters  we  have  considered  have  half  amplitudes  in  the  range  of  sorm  hundreds  of  kilometers.  This 
implies  that  the  contribution  due  to  the  imprecision  in  the  satellite’s  position  determination  (the 
ephemeris  error)  is  filtered  out.  Application  of  the  formulas  derived  in  [4] ,  yields  the  result  that 
ephemeris  noise  power  at  these  short  wavelengths  is  far  below  the  altimeter  noise  level  (~10  cm). 

We  summarize  below  some  of  the  uses  to  which  high  pass  altimeter  data  filtering  may  be  put: 

a.  Tectonic  detection: 

Seamount  detection 

Detection  of  signatures  of  elongated  features  extending  across  many  parallel  altimeter  data 
tracks.  These  signatures  may  be  so  weak  in  relation  to  the  noise  as  to  become  apparent  only  by 
correlation  across  tracks.  Examples  will  be  given. 

b.  High  pass  filtering  is  a  useful  tool  in  data  editing.  Isolated  data  errors  (perhaps  due  to  a  dropped 
bit  in  the  data  transmission)  and  of  small  magnitude  in  relation  to  the  unfiltered  data  (geoid 
height,  in  our  case)  are  made  to  stand  out  conspicuously.  An  automatic  data  editing  procedure 
might  ce  based  on  a  judiciously  chosen  filter.  Examples  of  data  error  detection  will  be  given. 

RECURSIVE  DIGITAL  FILTERS 

The  data  volume  to  be  processed  in  a  systematic  coverage  of  the  oceans  is  formidable,  and  therefore 
computational  efficiency  is  an  important  consideration.  Since  we  are  not  searching  for  a  well-defined 
signature,  there  is  no  well-defined  optimal  filter.  We  have  considered  the  digital  Butterworth  filters  [5 1 . 
They  have  the  advantage  of  great  simplicity,  and  their  frequency  response  functions  have  no  ripples. 
Applying  the  filter  recursion  in  a  forward  and  backward  sweep  removes  the  phase  [6,  p.  194) . 
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A  disadvantage  of  recursive  filters  lies  in  the  appearance  of  end  effects,  i.e.,  transients  near  the  ends 
of  the  data  span.  The  higher  the  filter  order  and  the  wider  the  frequency  pass  band,  the  greater  the 
end  effects.  In  Figure  1  we  have  plotted  the  frequency  response  functions  of  the  back  and  forth  sweep 
Butterworth  filters  for  some  typical  cases. 

The  Seasat  altimeter  data  are  available  to  us  as  sampled  values  spaced  at  about  3.32  km  (the  results 
of  averaging  over  500  separate  altimeter  readings). 

The  family  of  high  pass  filters  which  we  are  considering  involves  two  adjustable  parameters, 
the  order,  N,  and  half  power  wavelength,  X  (equal  to  the  half  amplitude  wavelength  for  the  back  and 
forth  sweep). 

The  proper  choice  of  these  parameters  depends  on  a  number  of  factors: 

a.  The  scale  of  the  tectonic  features  we  are  looking  for,  projected  along  the  subtrack.  For 
elongated  features  this  involves  the  angle  between  the  axis  of  the  feature  and  the  subtrack. 

b.  Control  of  data  gap  transients:  Due  to  the  finite  data  spans  (dry  land  and  other  causes), 
transients  are  excited  which  one  would  like  to  limit  to  a  small  neighborhood  of  the  data 
gaps.  For  fixed  X,  larger  transients  are  associated  with  high  order,  N. 

c.  Resolution  of  neighboring  signatures:  The  larger  the  order  N,  the  steeper  the  attenuation 
of  the  frequency  response  but  the  slower  the  decay  of  the  impulse  response.  This  leads 
to  poor  resolving  capability  for  high  order  filters. 

The  filtered  data  have  been  plotted  along  normals  to  the  subtracks  on  a  Mercator  grid.  This 
makes  these  charts  suitable  for  overlay  on  other  charts  (bathymetry,  basement  charts,  magnetic 
anomalies).  The  North-East  to  South-West  (descending)  tracks  have  positive  values  to  the  left,  the 
South-East  to  North-West  (ascending)  to  the  right.  The  normal  scale  has  been  chosen  so  that  1  meter 
of  geoid  height  corresponds  to  1°  of  longitude  in  the  grid. 

In  the  following  section  we  discuss  some  sample  results.  We  make  no  attempt  at  quantitative 
correlation  of  signatures  and  tectonics. 
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SOME  TECTONIC  SIGNATURES 

In  Figure  2  is  reproduced  a  bathymetric  chart  in  the  Mid-Pacific  Seamount  province  (7,  chart 
1803  N). 

In  Figure  3  one  observes  along  track  1 164  and  its  repeats  (thick  line),  a  “Seamount  Signature” 
at  12°  55'  N,  174°  8'  W  that  does  not  correlate  with  the  bathymetry.  Along  tracks  131  and  619  one 
notices  “Seamount  Signatures”  at  13°  55*  N,  175°  15*  W  and  14°  17' N.  175°  28' N  respectively, 
perhaps  due  to  the  same  feature,  but  not  correlated  with  the  bathymetry. 

A  weak  negative  signature  may  be  traced  across  several  tracks.  The  locations  and  magnitudes  are 
summarized  in  Table  1.  As  signature  magnitudes  we  have  taken  the  difference  between  the  central 
low  and  the  first  side  lobe.  Some  tracks  have  been  passed  through  the  wider  400-km  filter  frequency 
window  (see  Figure  1)  yielding  significantly  larger  signature  magnitudes.  Clearly,  the  200-km  filter  has 
not  captured  the  full  signature  and  we  have  a  rather  large  scale  feature.  However,  the  narrower  filter 
leads  to  a  better  separation  between  features  (resolution).  This  signature  may  be  associated  with  a 
slow  mesozoic  spreading  center  [8]  and  signify  a  trough  under  the  sediment. 


TABLE  1.  SUMMARY  OF  SIGNATURES  PERTAINING  TO  THE  LONG  FEATURE 

IN  FIGURE  6. 


Approximate  Low  Minus  High  (cm) 

Rev.  No.  Location 

N=2,X=200km  N=l,X=400km 


533 

11°41'N,  179°  26' W 

-40 

-70 

111 

11°  28' N,  179°  2' W 

-35 

-65 

619 

10°  20' N,  177°  6'  W 

-50 

131 

10°  9' N,  176°  48' W 

-30 

662 

9°  33' N,  175°  58' W 

-15 

-50 

1293 

9°  22' N,  175°  36' W 

-20 

There  is  also  a  weak  negative  signature  on  the  multiple  tracks  at  about  10°  32'  N.  1 75°  8'  W  of 
minimum  value  ~  -20  cm. 
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We  have  illustrated  an  important  point: 

Correlation  across  parallel  or  repeat  tracks  may  reveal  weak  signatures  not  otherwise  readily 
apparent  above  the  noise  level. 

In  Figures  4  and  5  we  show  a  bathymetry  chart  just  to  the  South  of  that  in  Figure  2(7]  and  data 
tracks  passed  through  the  N  =  3,  X  =  200  km  filter  (see  Figure  1)  respectively.  On  track791  at  4°  2l'N. 
172°  52'  W  a  “Seamount  Signature”  of  maximal  value  45  cm  is  noticeable.  It  does  not  correlate  with  the 
bathymetry.  The  large  excursions  at  the  data  gaps  on  tracks  102  and  260  are  due  to  filter  “End  Effects.” 
On  track  791  a  “Data  Error  Spike”  of -120  cm  is  noticeable  at  6°  32'  N,  171°  59'  W,  indicating  the 
utility  of  high  pass  filtering  for  purposes  of  data  editing.  This  is  essential  if  vertical  deflections  are  to  be 
estimated  by  automated  processing. 

In  Figures  6  and  7  we  have  bathymetry  [7]  and  filtered  data  tracks  (N  =  3,  X  =  200  km)  near  the 
Clipperton  Fracture  Zone  just  East  of  the  Christmas  Island  Ridge.  Along  the  multiple  track  (thick  line) 
there  is  a  “Seamount  Signature”  which  does  not  correlate  with  the  bathymetry.  Perhaps  the  ridge  to 
the  South-West  has  been  misplaced. 

We  also  notice  a  negative  signature  associated  with  the  Clipperton  Fracture  Zone  and  extending  over 
all  the  tracks.  It  may  indicate  a  trough  covered  by  sediment. 

Finally  in  Figure  8  we  show  some  tracks  through  the  Tasmanian  Sea.  Two  “Seamount  Signatures” 
are  apparent.  On  track  777  the  maximal  value  of  the  signature  is  56  cm  at  36°  19'  S,  159°  56'  E.  On 
track  533  the  maximal  value  is  49  cm  at  49°  8'  S',  1 57°  5'  E. 

CONCLUSION 

We  have  demonstrated  the  utility  of  high  pass  filtering  of  satellite  altimeter  data  and  presented  some 
results  in  chart  form  suitable  for  overlay  on  other  charts.  Full  interpretation  must  await  correlation 
with  other  data  classes,  where  available.  Satellite  altimetry  constitutes  the  most  extensive  data  class 
available  over  the  oceans  and  furnishes  important  constraints  on  crustal  models  and  bathymetry. 
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Figure  1. 


Figures  2,4,6. 


Figures  3,5, 7,8. 


FIGURH  C  APTIONS 


Frequency  response  of  some  high  pass  digital  Butterworth  filters. 


T  . Sampling  increment 

X . Wave  length 

to . Circular  frequency  =  2 ff/X 


Bathymetry  charts  taken  from  |7]  in  Mid-Pacific  Seamount  Province.  Mid-Pacific 
and  near  the  Western  Clipperton  Fracture  Zone  respectively.  These  charts  arc  to  he 
overlaid  by  Figures  3,5,  and  7.  respectively. 

Bands  of  subtracks  projected  on  a  Mercator  grid. 

Figure  3  for  overlay  on  Figure  2 
Figure  5  for  overlay  on  Figure  4 
Figure  7  for  overlay  on  Figure  6 

The  filtered  altimeter  data  have  been  plotted  perpendicular  to  the  subtracks,  so  that 
1  meter  corresponds  to  the  distance  between  adjacent  longitude  grid  lines  in  the 
Mercator  grid,  i.e.,  1  meter  ~  1°  longitude.  In  Figures  3.  5,  and  8  positive  values 
lie  to  the  left,  in  Figure  7  to  the  right. 
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Figure  2 
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